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Abstract 

We employ the influence functional technique to trace out the pho- 
tonic contribution from full quantum electrodynamics. The reduced 
density matrix propagator then is constructed for the electron field. 
We discuss the role of time-dependent renormalization in the prop- 
agator and focus our attention to the possibility of obtaining a dy- 
namically induced superselection rule. As a final application of our 
formalism we derive the master equation for the case of the field be- 
ing in an one-particle state in the non-relativistic regime and discuss 
whether electromagnetic vacuum fluctuations are not sufficient to pro- 
duce decoherence in a position basis. 
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1 Introduction 



Quantum field theory provides a broad framework witliin wliicli problems in 
diverse branches of physics can be formulated and addressed. Its success is 
primarily based on its inherent flexibility: its roots lie within both Hamilto- 
nian and statistical mechanics and this double aspect can be used invariably 
when dealing with particular physical phenomena. 

Its large domain of applicability in conjunction with the increasing inter- 
est in non-equilibrium phenomena has led in recent years into the adoption of 
successful techniques from standard non-equilibrium statistical mechanics or 
the quantum theory of open systems and their upgrading to fit the field the- 
oretic contexts. Hence the Feynman- Vernon influence functional technique, 
the Schwinger- Keldysh closed time-path formalism ( for extensive bibliog- 
raphy see reference |]l|) and Zwanzig's projection method p|, ^ have been 
used to deal with diverse issues including early universe cosmology, black 
hole statistical mechanics condensed matter physics and quantum optics. 

The earliest specimen among field theories is of course quantum electro- 
dynamics (QED). Being a theory of photons and electrons, it can be used 
to describe most of the phenomena in "ordinary" low energy matter. Hence, 
within its range of applications many issues connected with non-equilibrium 
phenomena arise. In this context one of the older (dating back to Lorentz 
for the case of classical electromagnetism) interesting questions is what the 
effective dynamics of charged particles when its interaction with its own 
electromagnetic field is taken into account. Or stating it differently, how the 
vacuum fluctuations of the electromagnetic field affect the evolution of an 
electron. 

More than that, the program of decoherence raises important questions 
relevant to the domain of QED. It has been proposed that the superselection 
rules for conserved charges could be another instant of the environment in- 
duced superselection rules |^ : that is, the electromagnetic field considered as 
an environment makes superpositions of states with different charge rapidly 
lose their coherence. Arguments in favor and against this proposal can be 
found in P, ||. 

Also of equal importance is whether single electron states exhibit deco- 
herence in position by virtue of their interaction with the photon vacuum 
[0, |[. The two-slit experiment suggests that this is not the case, but we 
would like to see in more detail the truth of this assertion. 
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These have been our physical motivations, but were not the only ones. 
We feel that the use of the influence functional technique in fleld theoretic sit- 
uations has been rather less frequent than the strength of the method would 
require. One reason for that might be the fact that it gives a description of 
the system in a Schroddinger picture, which has been rather rare in most field 
theoretic applications. Issues like renormalization (which is explicitly time- 
dependent) seem to be more complicated in this picture, Lorentz covariance 
is not manifest and calculations are definitely messier. 

Compensating for these problems arc two important benefits: 

1) Using the influence functional technique, we can more easily see the ef- 
fects of the initial condition of the environment in the evolution equations. In 
particular, the effects of the unphysical factorizing initial condition ( which 
is also used in the other techniques) are more easily understood and as we 
argue in section 3 we can find regimes where they can be removed by a renor- 
malization process. 

2) The main object in the influence functional formalism is the reduced den- 
sity matrix propagator, the knowledge of which enables us to construct mas- 
ter equations, study the classical limit, decoherence phenomena and pass 
from a field description to a particle one. 

From the above discussion, we believe that our aim in this paper has been 
made clear. We intend to use the influence functional technique to trace 
out the effects of the electromagnetic field in full QED, hence obtaining 
the reduced density matrix propagator. This is meant to be our tool for 
examining the possibility of a dynamical origin of the charge superselection 
rules. As a special application we shall construct from this the evolution 
equation for the regime where the initial field state lies in the one - particle 
sector, aiming to obtain a master equation for a single electron as it interacts 
with the vacuum fiuctuations of the photon field. 

In the next section we give a brief review on the infiuence functional 
formalism, fixing primarily the notation and introducing some important ob- 
jects. The calculation of the propagator for the free spinor field with external 
sources (of great importance for perturbation theory) is also briefiy sketched. 
Section 3 is the largest and from the technical side the most important one. 
Herein the reduced density matrix propagator for the electron field is cal- 
culated, using perturbation theory. The section includes a discussion of the 
factorizing initial condition, how one can remove the inevitable divergences 
through a renormalization procedure and the implications of our result for 
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dynamically induced superselection rules. In section 4 the special case of a 
single non-relativistic particle is studied. A master equation is obtained and 
the question of whether the photon vacuum can play the role of a decohering 
environment is studied in detail. 

Before proceeding any further we should mention that a master equation 
in the context of quantum electrdynamics for the reduced dynamics of the 
charges has been derived before It is, to the best of our knowledge, the 
first attempt to derive such an equation in QED but it does not utilize the 
functional Schrodinger picture and it focuses on different issues. 

2 The model 

Our starting point is the QED Lagrangian 

i){%^^dm - m)i) - \f^uF^^ + eij^A^i; (2.1) 

We are going to use the influence functional technique to trace out the photon 
field and derive the master equation for the reduced dynamics of the electron 
field. To do this we shall work in the Bargmann representation (considering 
essentially coherent state path integrals) and under the Schroddinger pic- 
ture. For more details in the various elements we employ in our formalism 



the reader is referred to |T^, [IT], 0. The relevant object is the evolution op- 
erator e~*^*, the matrix elements of which can be read from its path- integral 
expression 



(^/, A^le-^^^IV'., A>) = J [#][#] [rfAle^^^'^'^^'^l (2.2) 

where the sum is over fields such that ip{0) = ipi, iplt) = tpf, A^^lO) = Af 
and A^^lt) = A'^. The knowledge of this object is sufficient to determine the 
reduced density matrix at time t 

Pit) = Trph (e-^^Voe*''*) (2.3) 

where the trace is over the photon field degrees of freedom. The reduced 
density matrix propagator J is defined by 

Pi(V'/, ^'f) = I D^lD^.Ji^f, i,'- t\i,[, 0)po(^:, ^i) (2.4) 
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and can be written in a path integral form 



\dij'][d^']exp \ iSS,ij] - tSeiij',^'] \W[^,^lj;ij',ij'] (2.5) 



with integration over paths ijj{t) = ipf, ip{0) = ipi, V''(0) = ipi and ip'{t) = ip'f- 
Here, W is the influence functional containing the information of the 
dynamics and initial state of the photon field. In the case of factorisable 
initial conditions po = Peo ® PpW this has a path integral expression 



W[i,,^-i,\^'] = J DAfDADA[ J[dA][dA']exp!^iSph[A]-iSph[A'] 

exp\iSint[A,ip,'^]-iSint[A',ip','^'] \ x ppho{Ai,A'i) (2.6) 



with path integration over vector fields such that A{0) = Ai, A'{0) = A'- and 
A{t) = A'{t) = Af. Since the Lagrangian contains an interaction term of 
third order to the fields, we have to rely on perturbation theory. For this it 
is sufficient that we compute the spinor propagator with external sources. 

2.1 The spinor propagator with external sources 

Let us consider the Lagrangian density 

C = '^p(i'y^'^'' - m)^ + f]^) + iprj (2.7) 

in terms of the Grassmann fields ip, ijj and sources t], fj. Written in terms of 
the Fourier transformed fields and sources 



^ = Ei/pe-''' ^^ = T.yy'^ (2-8) 
p p 



p 



^ We denote the path-integration measure (over spacetime fields) by [d.] and the mea- 
sure over fields at one Cauchy surface (single moment of time ) as Z?; the former is a 
formal expression while the latter is the Gaussian measure with respect to which the field 
Hilbert space is defined. 
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it reads 

L = / d'xC = yli^^o - il'YPi - m^')y, + r)lj\ + yl^\ (2.10) 
which corresponds to a Hamiltonian 

^ = E yli^'fPi + ^7°)2/p - vliX - yh\ (2.11) 
p 

Now, Up is a four-component Grassmann variable. In order to write the 
Hamiltonian as a sum of single-component Grassmann variables one has to 
diagonalise the matrix 

7°7>i + m7° (2.12) 

Writing in the standard way (suppress the index p) y = yiUi + ziVi, rj — 
TjiUi + r/'wj, i = 1,2, there is a particular choice of u^v such that the above 
matrix becomes diagonal. The Hamiltonian becomes therefore a sum of terms 
of the form 

H = ujyy -fjy-yr] (2.13) 

where uj^ = (p^-|-m^). This is just the Hamiltonian of a Grassmann harmonic 
oscillator with external sources, the propagator of which reads in a path- 
integral form 

(y/,i|y,,0) = J DyDye'^^y'y^ (2.14) 
where summation is over paths with y{0) = y^, y{t) = yf and the action is 

iS[y, y] = yy{t) + i f ds{iyy - H{y, y)) (2.15) 
Jo 

The path integral integral is readily evaluated by the saddle-point method 
to yield 

{yf,t\y,,Q) = e'^^^ (2.16) 

with 

iSd = yfyic-"^' + I f ds[e-'^^'-%rj(s) + e-'^%s)yi] 
Jo 

-[ ds r ds'e-'"^^'-''^ r]{s)7]{s) (2.17) 
Jo Jo 
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One can thus use the inverse Fourier transform to construct the propagator 
for the full spinor field 



{ipf, t\ipi, 0)nfj = exp <^ I d xd X 
+i 



ipf{x)A{x,x';t)ipi{x') 



I ds(^^f{x)A{x, x'\ t — s)r]{x', s) + r]{x, s) A(x, x'; s)4'i{x') 

t i-s' 

ds I f]{x,s)A{x,x']\s — s\)'r]{x ,t) 
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where A is a "propagating kernel" 



(2.18) 



A{x,x'-t) = I ^e-^P("-"') (" ^^7° ^>- cosu;„t-2sinu;„tl (2.19) 
J [2'kY \ m J 





3 The reduced density matrix propagator 

Let us concentrate now on the photon field dynamics. The relevant part of 
the Lagrangian density is 

C = --^F,,F^'' + J,A>' (3.1) 

where = ip'^^ip 

The main idea in the influence functional formalism is to proceed and use 
the propagator for the above Lagrangian , while treating J as an external; 
source. Now this is clearly a linear action, hence the propagator in presence 
of sources can be computed exphcitly. 

It is more convenient to fix the gauge for the photon field {d^A^ = 0, 
= and expand the spatial components of the field in normal modes 

Ax) = E E {qM coskx + g^,(k) sin kx) (3.2) 

k r 



is essentially a unitary operator in the single-particle Hilbert space out of which 
the field Fock space is constructed. 
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where the polarization matrices satisfy 

^e^(r)(k)e^(r)(k) = T'^{k) = 5'^ - k'^/k"^ 



(3.3) 



Substituting this into (2.6) we get the following expression for the relevant 
part of the action 

SpH + Si^t = I f dsJ2Y.Y.(€i(^) + k'^'M + 2gw(k)Jw(k) (3.4) 

^ k r I 

In this expression the index / takes values 1 and 2 according to the decom- 
position (3.2) and the external fields J read 



Jij.(k) = e y o?^xcoskxJ*(a;)ei(r)(k) 
J2r(k) = e y d^xsinkxJ*(a;)ei(r)(k) 



(3.5) 
(3.6) 



The action corresponds to an infinite number of harmonic oscillators un- 
der external source and the propagator can be computed the standard way. 
Moreover if we assume the initial state factorisable and the photon part to 
correspond to the vacuum = |0)(0| we can use equation (2.19) to compute 
the influence functional. It is a standard calculation to show that 

W[J, f] = C{t)Y[exp \-i fdSi r dS2[Ja - J'alisi) 

L Jo Jo 



X 



2ujk 

2 



sina;fe(si - S2) [Ja + J'a\is2) 



(3.7) 



where a is a collective index containing k, I and r. Using the inverse Fourier 
transform we obtain the final result for the influence functional 

W[J, J'] = C{t) exp j d^xd^x' dsi J^' ds2 
- J"]{x,s,)riij{s^ - S2,X,X)[J^ + J'^{x',S2) 
-[J' - Ji]{x,Si)l^ij{Si - S2,X,x)[Jj - J'j]{x',S2))] (3.8) 
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Here, rj and u are the dissipation and noise kernel respectively which read 



V^As,x,x) = e'J ^-fL_e-^^^^-^^r^{k) sin 00,3 (3.9) 

iy,j{s,x,x) = e^ [ -^^e"*^(^-^)r^(k) coso^^s (3.10) 

Having the expression for the influence functional we can proceed to eval- 
uate the reduced density matrix propagator. But first we should recall that 
its derivation depended upon the assumption of factorised initial state. We 
should first then discuss what is the physical meaning of such a condition. 

3.1 The initial condition 

Our choice of the initial condition (factorisable density matrix, vacuum for 
the photon field) has been such as to make the calculations easier, but on 
the other hand it introduces a number of problems . This condition has 
been extensively used in quantum Brownian motion models and is generally 
deemed unphysical. Such a separation between system and subsystem re- 
quires a large amount of energy (in the case of fields an infinite amount) and 
hence cannot be taken as realistic. The effects of this initial condition have 
not been fully understood yet. It was conjectured [1^ that it is related to 
a jolt in the diffusion coefficients at early times, dissapearing in a timescale 
proportional to A~^, where A is the high frequency cut-off of the environ- 



ment. Nevertheless this is not probably the case. Recent studies |TJ] have 
shown that this jolt is present even if one assumes correlated initial states. 
The typical behavior in those systems is that the preferred degrees of freedom 
at early times couple very fast to the high frequencies of the environment, 
until a dynamical equilibrium is established and then at later times the low 
frequency modes start playing an important effect (in essence the long time 
limit is governed solely by the latter) . In this class of systems it is usually the 
case that uj where uj is the typical frequency of the preferred degrees of 
freedom. To conclude, it is the existence of the above dynamical equilibrium 
allowing us to trust the analysis with the given initial condition. 

Doing field theory just makes things more difficult as one would expect. 
In QBM models the separation between system and environment follows 
naturally by assuming the typical timescales of the problem and somehow 
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the environment appears to be "robust" to the relevant degrees of freedom. 
Backreaction effects need not be taken into account as far as the environment 
is concerned. One of the concequencies is that a cut-off naturally appears 
only in the frequencies related to the environmental degrees of freedom (A). 
In QED this separation is not obvious. We could equally well formulate the 
inverse problem i.e. what happens to the photon field if we trace out the 
fermionic degrees of freedom? This, for example, is a question of relevance to 
solid state physics. Still we expect that some of the features of QBM models 
will be relevant but we have to be very careful in our analysis. 

Our aim is to find the regime where our initial condition can be a good 
approximation, in the sense that a suitable renormalization of the parameters 
can be expected to give reliable results. What we assume is that at times 
t < the electron - photon field combined system lies on its ground state 
and at time t = an operation is performed on the spinor sector, which 
effects to a local excitation of spinor degrees of freedom. Our aim is to see 
how the reduced density matrix describing these degrees of freedom evolves 
in time. The vacuum of the full QED theory is a correlated one, so no 
operation carrying finite energy is sufficient to fully separate the spinor from 
the electromagnetic degrees of freedom. There will be correlations stemming 
from the ultra- violet and infra-red sector of the photon spectrum. 

Now, in perturbative quantum field theory, one is starting from the 
Hilbert space of free fields (hence considering a factorised vacuum state) and 
then the renormalization procedure provides a mapping (though a mathemati 
cally ill-defined one from the free fields Hilbert space to the interacting 



one ITSf. This observation suggests a rephrasing of our problem : is there 
any renormalization procedure enabling us to remove the effects of the initial 
condition? 

The expected behavior is the following. At early times the high frequency 
modes of the photon field couple very fast to the slow ones of the fermion field 
and vice- versa. This latter coupling expresses the backreaction of the relevant 
degrees of freedom to the environment (We believe backreaction effects to 
be important as far as decoherence is concerned). Our experience tells as 
that it will only affect the initial moments of the dynamics of the photon 
field had we wished to study them. Hence in the times where the problems 



Since the free and interacting field Hilbert space carry unitarily inequivalent repre- 
sentations of the canonical (anti)commutation relations 
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of the unphysical initial condition are expected to arise the contribution of 
correlations with soft photon modes is not important. This leads us to expect 
that the standard renormalization procedure of the propagator must be 



sufficient to account for the effect of the soft modes. This will be verified 
in our analysis, that no divergent contributions stemming from low energy 
modes appear in the density matrix propagator apart from the standard ones 
due to perturbative evaluation. 

The case is more complicated for the highly energetic modes, since there 
exists an infinite energy barrier between our factorised initial state and a 
physical one {i.e. the Hilbert space of the interacting theory does not contain 
a factorized state). This implies that the time it will take for these modes to 
come to equilibrium will tend to become infinite. To get over this difficulty we 
may observe that no such modes can be excited with energy larger than the 
total energy of the system. Hence there exists of a natural cut-off A (state- 
dependent) for both spinor and photon modes (This is due to the fact that 
photon and fermion fields involve typical frequencies which are of comparable 
order of magnitude. There is not a preferred time scale distinguishing one 
from the other in a natural way). Introducing this cut-off can be taken as 
considering factorised initial state for only the modes with energy less than A, 
while higher energetic ones (for both fields) can be taken to be in dynamical 
equilibrium This implies that the effect of the spurious initial condition 
is to be contained in the parts of the propagator containing the cut-off. 
Hence a renormalization procedure (additional to the standard dynamical 
one) removing these terms would yield reliable results for sufficiently large 
times. 

More formally, one could write that the physical density matrix should 

be 

P0 = Pe® Pph + e^PcoT (3.11) 

The first term contains a factorized state for all modes it describes while 
the second one gives the contribution of the correlated vacuum of the higher 
modes. In other words we are asking for a state which in the low energy 
modes is essentially uncorrelated while in the high energy ones contains the 
correlations of the real vacuum state of the interacting theory. The later 
contains differences from the factorized vacuum of the order of e^. Time 



^ More precisely we can assume a truncation of the field theory such that high energy 
modes do not contribute in the dynamics and hence lie close to their vacuum state 
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evolution to lowest order in perturbation theory takes the low energy modes 
to low ones and similarly for the high ones. To have mixing of the regimes 
we must consider orders . This means that the contribution of the initially 
correlated terms (already of order e^) to the low energy ones at late times will 
be of order e^, hence we can safely ignore it in the lowest order of perturbation 
expansion. 

To summarize: there is a regime of energies E, where a (time-dependent) 
renormalization on the propagator is expected to remove all effects due to 
our unphysical initial condition and our expressions can be valid for times 
t » E~^. In particular, the non-relativistic limit for a single particle where 
E c:^ m should definitely lie within our range of validity. 



3.2 The propagator 

The index notation It is more convenient to use an index notation for 
the spinor fields so that we will not have to carry the explicit dependence on 
X. We have essentially four types of spinor fields: ip{x), ip{x) , ip'^x) ,ip'{x). 
The first two propagate forward in time and the other backwards (note the 
similarity with the CTP formalism). Now, ip and ip' in the path integration 
are fixed by their initial condition at t = 0; we shall assume that they carry 
an upper index (in place of their x dependence), while ip and ip' are fixed 
by their values at t and we shall use a lower index for them. This is an 
elaboration of the notation used in [Q. 

Contraction of indices amounts to integration with respect to x, while 
kernels carry indices according to the fields they are contracted. Hermitian 
conjugation for a kernel amounts to inversion of all indices. Hence for exam- 
ple our propagator (2.19) reads in this notation 

{ijf, t\ij,, 0) = exp \ijfaA{t)tij'i + I f dsiijfaAtit - s)r]\s) + 





Us)At{s)^',) - fds f ds'Us)Al{\s-s'\)^\s') 
Jo JO 



(3.12) 



Whenever an upper and lower index do not denote integration and are to be 
thought as free a bar will appear in top of them. 



Perturbation expansion To compute the propagator one needs to eval- 
uate the path integral (3.8) using the expression (2.6) for the infiuence func- 
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tional. Since the influence functional contains term of fourth order to the 
fields, an exact calculation is infeasible, hence we will have to rely on pertur- 
bation theory. For this we shall employ the following identity 



{i:f,t\iJi,o)^^{ij'f,t\i:'-,o)^,-, 



\ri=ri=ri' =ri' =0 



(3.13) 



Expanding W in powers of one can use this formula to derive the pertur- 
bative series for the density matrix propagator. 

Now, in the expansion we shall encounter loop terms (corresponding to 
Tr A(0) or A"b(t — s)A^a(s). Their effect would be to produce a change on 
the kernel A and to necessitate a field and coupling constant renormalization 
. We will prefer to work in the zeroth loop order, that is consider only tree 
diagrams [| . 

Hence, when taking only the tree diagrams into consideration we see that 
in the propagator there appear exponentiated four interacting terms. 



— I ds I ds' 



+ ds ds' 



+ ds ds' 



— ds ds' 



i:fcA'-g{t-s')jjA^,{s')ij. 
^faA\{t - s)YAM;^{s)f,:^^,^{s - s') 

i,':Aj^{s')^^A/{t-s'Wf, 
i,'':A,\s)fA,\t - s)fj,f^,,,,g{s - s') 

^fA',{t - s')7^A5,(s')^f 
^P':Aj{s)YA,\t - 5)^>*,,e,(5 - ^0 
^'^A/(.')7^A/(t-/)^;, 



(3.14) 



^We could obtain the same result had we performed a saddle point approximation in 
the path integral and used perturbation theory in order to construct the corresponding 
"classical" equations of motion. 
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In the above expression fi{s — s') stands for {u + iri){s — s') and fi* denotes 
its complex conjugate. We observe in these terms combinations of the form 
A7A, which amount to products of three and two 7 matrices. For physical 
interpretation it would be more convenient to write these in terms of a fixed 
basis of the 4x4 vector spaces where the Grassmann variables. Actually in 
our case (parity preservation) only 7^*, 1 and a'^'^ are relevant. One can use 
the identities 

^M^^ = r]'"'l + m^"' (3.15) 
= g^'Y + 9"'^ - (3.16) 

to bring all interaction terms into sums of the form {^jjAilj){iljBip) where A 
and B are matrices: 1,7'*, a^". 

Let us denote by capital indices I.J, . . .whether the indices correspond 
to scalar S, vector V, tensor T current (in terms of Lorentz indices S stands 
for no index , V for a single index and T for two antisymmetric indices). 
With this notation we can write the reduced density matrix propagator as 

j(Vi/,V'/;^IV'i,V^;,o) = 

e^.a'/^'iAi^'f.i^'.Aji^'f, - Cu^c^aAii^'fi^'lAji^'f,^ (3.17) 

The coefficients e and ( can be exphcitly computed. Their Fourier transforms 
are 

ejj{x,x',y,y') = / ^i!^e^p(--')-P'(.-.'),^^(p,y, (3.18) 
Cjj{x, x\ y, y') = I ^|!^e^^'(--')-P'(.-.')^^_^(^,y, t) (3.19) 

Renormalization Divergences appear into the coefficients of the propaga- 
tor. We have chosen a naive reqularisation scheme of substituting / d^k with 
/ (Pke~'^^. Divergences then take the form of either poles in e or logarithms 
of e. The renormalization procedure we are going to follow consists in essen- 
tially dropping all terms containing the divergences. There is an ambiguity 
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as far as the logarithmic divergences go, in the sense that one can with a sim- 
ple rescaling of e terms with no dependence on the cut-off appear. Happily, 
in the case of such divergences the ambiguous terms are oscillatory in time 
and generically negligible at late times, when the approximation is thought 
to be valid. 



3.3 The master equation 

From the form (3.17) of the reduced density matrix propagator it is straight- 
forward to obtain the master equation. Recall that the relevant procedure 
consists in writing all multiplicational action of initial variables operators on 
J to differential action of final variables. It can be easily shown that the 
master equation is of the form 

r d^p d^p' d^k f . ^ , , , ^ , , x ? / / , \ -^n 
+ 1 (^(^(^(^^e "/./(P'P'k;^)['^/(P>k)Jj(p,k),p] 

-e2/?,j(p, p', k; t)[J,(p, k), [Jj(p', k), p]]) (3.20) 

The kernels a and /3 can be computed from the knowledge of e and 

It is important to remark that the interaction of the spinor field with 
its environment comes through terms quadratic in the currents J/(p,k) = 
^Ai^{p^k). 

What are the consequences of the expression (3.20.) as far as decoherence 
phenomena for the electron field are concerned? The complicated expressions 
for our coefficients does not support the probability of giving a general ar- 
gument or proof. What they hint is that if we are to deal with those issues 
we must restrict ourselves in particular regimes. Such a regime (of a single 
particle state in the non-relativistic limit) we shall study in the next section. 

Still, based on the experience obtained by study of simpler open systems 
one can make a number of conjectures with high plausibility. Both 
dissipation and diffusion terms in the master equation are expressed in terms 
of the three currents (scalar, vector, antisymmetric two-tensor). This means 
that they form the channels through which the spinor field interacts with 
its environment. In particular, the diffusion term for each current suggests 
that the reduced density matrix has the tendency to become approximately 
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diagonal in some basis close to eigenstates of the current operators As it 
stands, we cannot determine which such a basis can be, let alone whether the 
suppression of the off-diagonal terms is sufficiently strong. But if such regimes 
exist, then superposition of states with large difference in the expectation 
values of the current operators would decohere fastly. We could then say that 
we have some sort of "local current " environment-induced superselection 
rules. This is a more general case than a single charge superselection rule and 
actually seems more physical: Charge superselection rule is a special case of 
this, but we could also have decoherence between states with the same global 
charge but different local charges, i.e between two seventeen-electron states 
with expectation values of currents having support in spatial neighborhoods 
with large separation. The study of such regimes would be invaluable towards 
understanding the role of decoherence in quantum field theories and can in 
principle be performed starting from the propagator (3.17). In connection 
to the above it is also worth noticing that not all of the currents appearing 
in Eqn (3.20) are conserved. Only the vector current has got this property. 
This current is thus related to a conserved quantity, a charge, given as a 
volume integral over the corresponding local density. When this volume is 
taken to infinity it will be an exactly conserved quantity. Then in this case 
according to we have exact decoherence which because of the conservation 
law is related to the existence of an exact superselection rule. Nothing like 
that is obvious from the form of the above master equation. It is not very 
clear to us how the results of ^ could be reproduced in the given formalism. 
We believe it to be an issue worth investigating. Another interesting point 
could be made in relation to the ideas of Gell-Mann and Hartle that in 
decoherent histories hydrodynamic variables are the variables that habitually 
decohere, being the quasiclassical variables of genuinly closed systems. These 
variables are related to locally conserved currents. Since her not all of our 
currents are conserved there is an open question of what decoherence in the 
basis of those currents (if the phenomenon occurs) means and whether the 
basis of the conserved current is in any sence prefered over the other ones as 
one would naively expect. 
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4 The master equation for a single particle 



The reduced density matrix propagator for the spinor field (3.17) is quite 
general, and can act as a starting point for exploring particular problems . 
In this section, we will derive the master equation for a single non-relativistic 
particle, having as our starting point the field propagator. 

The general technique for the derivation of master equations for particle 
dynamics, starting from field theoretic evolution equations has been devel- 
oped in reference 0, the notation of which we shall employ. In our treatment 
we are going to make the following assumptions: 

i) The initial state of the field lies within the one-particle sector. 

ii) We work in frame of reference such that the momentum of the initial 
state lies in the non-relativistic regime: Ip/^nl << 1. This means in particu- 
lar that there is no pair creation and the master equation is expected to be 
trace-preserving. 

iii) We work to the lowest non-trivial order in and |p/'ni|. 

The assumption (i) guarantees that the propagator of the single particle 
is of the form PJP, where J is the field propagator and P projector on the 
space of field states that project to the one-particle sector. It reads explicitly 

Pp = PipPi (4.1) 

where the operator Pi is defined by its symbol in the Bargmann representa- 
tion 

Pi{^,^)=^P^ (4.2) 

Here P is a 4 x 4 matrix that projects into the two dimensional subspace of 
Dirac spinors that correspond to particles. In the Dirac- Fermi representation 
for the 7 matrices it reads P = diag{l, 1, 0, 0). 

The propagator Jiparfor the single particle then reads 

x(x/P^;)(^;Px})J(^/.^};t|^:,^i;0)(x:P^:)(V5.Pxi) (4.3) 
It is easily shown that Jipar is of the form 

JipariXf, Xf, t\x\, Xi-, 0) = K\/xfa X/ x[" Xid (4.4) 



17 



The kernel K^'b/ = Kr.r sfs- {xf,yf;t\xi,yi) is the single particle propagator 
( the indices r and s take values 1 and 2 and correspond to the two spin 
polarizations). To lowest order in the kernel reads 



be 



+ [ ds [ ds' 
Jo Jo 



+ 



f ds r ds' 

Jo Jo 



PA'^.it - s)yA%(.)P/.^ - s'] 
PAr{s')^^Aj{t - s')P 

PA\{s)YAMt - s)Pflij,eg{s - S')] 

PAj>{t - s'HA/{s')P 



(4.5) 



The non-relativistic limit is more clearly obtained in the Dirac-Fermi repre- 
sentation for the 7 matrices 



7 



1 
-1 



a' 
a' 



In this representation the kernel A reads (its Fourier transform) 



A(p;t) 



— cos (jJrit — i sin ujrit — — sin u!„t 

— —smcujjt —— cos curit — i sin u!j)t 

ml' m f y 



(4.6) 



(4.7) 



Substituting this into (4.5) we derive the following expression for the propa- 
gator 

K (Xf Vftlx- V ) = I ^""P ^^^P' ^i{p\xf-Xi)-p{yf-yi)-^{p^-p'^)) 



(27r)3 (27r) 

X 1 + / -^e'^^'^^-y^^T^ik) t ds r ds' 
w? J (27r)3 ^ ' Jo Jo 

xA'it, s){a''a%^r,e-''^^-''\a^a%^sA^''{t, s')y 
x{A'{t,s)r{a'a%^s.e''^'-^'\a^a%^sA'\t,s') 



(4.8) 



where the functions A denote, 

A'it,s) = tp'f{t,s) + k'git,s) 
A''{t,s') = t/f{t,s') + k'g{t,s') 



(4.9) 
(4.10) 
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with, 



f{t,s) = e-'"^' + cos m{t- 2s) 

9(^,8) 



^-im{t-s) g-^^^ 



(4.11) 

(4.12) 



and ( )* means complex conjugation. Simplifications occur, when we ignore 
the spin degrees of freedom, that is we consider initial and final states of 
total ignorance for the spin. Hence we define the "spinless" propagator 



K[Xf,yf,t\x„y,) J (27r)3(27r)3'' 



)) 



X 



1 + 



m? J 



^ik{xi-yi)rpij 



) ds ds' 
Jo Jo 



(27r)3 



x(A^(t,s))*e^'=(-^')(aV%,,A''(t,5') 



We now perform the k integration. Writing r = Xj — we get 



(4.13) 



d^k 1 
(27r)3 2A?^ 



-ikr—iks 



1 



47r^ — 



Svrr 



6{s-r)^I{r,s) (4.14) 



In this integral the principal value is assumed for the real part, and a delta 
function on advanced time s + r has been dropped out since we have positive 
s in our integration range. 
Similarly we can obtain 



/ 
/ 



(27r)3 2A; dn ^ ' ^ 



d^k 1 
(27r)3 2fc' 



-ikr—iksj^2 



Qj.2 



I{r,s) 



(4.15) 
(4.16) 



Using these expressions it is straightforward to perform the p and p integra- 
tion to obtain our final expression for the propagator 



K{xf, yf, t\xi, yi) = C{t) exp 
W[a{t,r){xf - Xi){yf - yi) 



im 



"-^[{xf-Xif - [yf-yif] 
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if3{t,r) 
H — [Xf - Xi){Xi - Vi) 

-i(3*{t,r){yf - yi){xi - Vi) -'y{t,r)] 



(4.17) 



The expressions for the coefficients a, (3, 7 are quite comphcated. What we 
are more interested is their dependence on r. RecaUing that we are working 
in the non-relativistic regime and for the dimensions of our system r can be 
taken much smaller than t. Hence to a first approximation one can expand 
the real part of /(r, s) in a power series with respect to r , which essentially 
corresponds to a power series with respect to 1/c (not so with the imaginary 
part due to the presence of the delta function) . This gives for the coefficients 
following r dependence 

a{t,r) = ai{t,r)/r + A{t) + 0{r) (4.18) 
(3{t,r) = ^bi{t,r) + 0{r) (4.19) 

= ^7i(t,^) + 0(r) (4.20) 

where 71 denotes the integrals containing the delta function which are 

calculable in terms of elementary trigonometric functions (see appendix) . 



4.1 The Gaussian (Coulomb ) approximation 

Still, the r dependence of the propagator is very complicated. Things sim- 
plify if one considers the "Coulomb" limit c — > cxd or r — > 0. This essentially 
corresponds to ignoring the retarded propagation effects for the virtual pho- 
tons and is expected to be a good approximation (at least for early times) 
when the spread of the wave packet is very small. Ignoring for the moment 
the domain of validity of this result, we can easily see that the functions f3 
and 7 vanish and a{t) contains only a time dependent contribution., hence 
the propagator becomes Gaussian in this limit 

K{xf: yf t\xi: Vi) = C{t) exp 

+e^{a{t){xf - Xi){yf - yi) 



(4.21) 
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From this with the standard procedure one can construct the master equation 
for the reduced density matrix of the "spinless" particle: 



d_ 
di 





(4.22) 



or in operator form 




d 



2e2 



(4.23) 



m 



where 



m 



(4.24) 



and its plot with time is given at the figure 1. We can see that it is oscillatory 
with a basic frequency of the order of m~^. It is rather straightforward having 
this simple expression for the propagator, to examine whether in this regime 
there is any decoherence in position for the electron. For this we consider 
the evolution of the initial state 



which is a superposition of two wavcpackcts, with their centers separated by 
L. Time evolution with the propagator (4.12) yields a mixed state with two 
diagonal terms propagating according to the classical equations of motion 
plus two off-diagonal terms. To examine the appearance of decoherence one 
is interested in computing the prefactor of the latter. This is found to be 



It is easy to check that the off-diagonal terms are not suppressed. A necessary 
condition for that would be a growth of f3(t) of at least t"^ with time, while 
we can see from figure 2, that j3{t)/t^ is falling asymptotically as 1/t. So 
within this approximation , no decoherence in position seems to be possible 




(4.25) 




(4.26) 
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Figure 1: The function f{t), containing the effects of diffusion at times t » 



0.075 
0.05 
0.025 



-0.025 
-0.05 
-0 . 075 




22 



Figure 2: 



The function I3{t)/t^ for t » m-\ 
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for the electrons. This is in accordance with the facts of the classic two-slit 
experiment. 

But how reliable is the approximation of switching off the r dependence 
of the propagator, or rather in which regime should it be trusted. There are 
two remarks we should make. First, the term containing 7(t,r) is increasing 
with time, hence at later times (depending primarily on the initial configu- 
ration of the system) this term becomes sufficiently large , to invalidate our 
perturbation scheme. On the other hand in the Gaussian approximation the 
perturbative terms in the propagator remain sufficiently small at all times. 
We therefore conclude, that as long as our Gaussian approximation is valid, 
perturbation theory is reliable. Of course, both the Gaussian approximation 
and perturbation theory is due to break at the long time limit. 

What then can we say, about decoherence in position of the electron, 
outside this approximation. Unfortunately very little, because the term that 
we would expect to give rise to decoherence (namely 7(t,r) ) is the one 
that eventually breaks perturbation theory and hence docs not allow us to 
reach any general conclusions. The problem is aggravated, because the only 
natural length scale in our system is the Compton wavelength of the particle 
and qualitative arguments of length scale separation are not sufficient to give 
any insight. 

What we can conclude from our analysis is the following. There exists 
a regime dependent on the spread of the initial state of the electron, that 
one can safely use the Gaussian approximation. In such a regime no deco- 
herence phenomena appear. As the spread of the initial state L increases, 
the influence functional phase gets a large real negative contribution which is 
characteristic of decoherence phenomena, and we are led to conjecture that 
decoherence is indeed possible. 

One can view these results in the light of the following qualitative ob- 
servations. Decoherence in any system is connected with the propagation 
of coherent phase from the system to the environment. In our particular 
case, the carrier of the information is of course the photonic field. Hence 
decoherence is expected, when the configuration of the state of the particle 
is suitable to excite a number of photons , which propagate in the environ- 
ment carrying the phase information. In the Gaussian regime as we have 
seen the flow of photons (essentially dissipation) is absent, while the (3 term 
in the propagator- negligible in this case- essentially contains the effects of 
dissipation. 
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Now in the electron - photon vacuum system there are essentially two 
competing processes governing the emission of photons. On one hand the 
space where the electron wave function has support "contains" a charge dis- 
tribution, which essentially acts as "screening" medium, that does not allow 
virtual photons to escape. The vacuum fluctuations here are of an essentially 
" Coulomb" type and can only give noise, whose effect only effect is random- 
ness in the evolution of the electron. This process is stronger when the size of 
the wavepacket is small, and it is precisely in this regime that our Gaussian 
approximation is valid. 

As the spread of the wavepacket increases, and in particular its shape 
deviates more from a spherical distribution another process starts gaining 
momentum. Higher moments of the charge distribution are coupled to the 
electromagnetic vacuum and as such they enhance photon emission. In par- 
ticular, the case where the electron is in a superposition of two spatially 
localized states, the contribution from the dipole moment (rather than the 
quadrapole for more localized states) becomes significant and the electron 
essentially behaves as an oscillating dipole. When this effect becomes strong 
and sufficient to overcome screening, significant photon production occurs 
and decoherence is possible. 

At this point we should comment on the relation of our results with the 
existing bibliography. The density matrix propagator is distinctly different 
from the one obtained in [^. These authors have been able to derive a 
Markovian limit of their master equation and obtain the Lorentz equation in 
the semiclassical limit. The crucial difference is that they have derived their 
result under the dipole approximation. For this they considered a restor- 
ing harmonic potential so that the particle is constrained to lie in a small 
spatial neighborhood. In our case having derived our master equation from 
full Quantum Electrodynamics, we have considered an essentially free elec- 
tron. Since the only natural length scale appearing in our calculations is the 
Compton wavelength of the particle, and free evolution does not allow for 
localization of the particle, the dipole approximation should not be expected 
to be valid here. 
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5 Conclusions 



The density matrix propagator (3.17) and (4.19) for the single particle , we 
consider to be the most important results in our paper. The former is to 
be seen as a starting point for many interesting applications (to be outlined 
below), while the latter provides (we believe) a conclusive demonstration 
that electromagnetic field vacuum fluctuations are not sufficient to cause 
decoherence to single free electrons. 

It has been rather disappointing for us not to be yet able to say anything 
conclusive yet about the possibility of dynamically induced superselection 
rules. As we said earlier the complicated nature of our expressions preclude 
any possibility of being able to make any statement about generic states. Still 
we think we have provided the basic tools for a more thorough investigation 
of these issues. 

The study of various regimes seems to be a natural step after this work. 
An important case can be the case where the field lies either in the or 
the M particle state. The propagator would then reads PnmJPnm where 
Pnm is defined by PnmP = {Pn + Pm)p{Pn + Pm)-, Pn being the projector 
in the A^ particle sector. If the hypothesis of the dynamical origin of charge 
superselection rule is valid then decoherence phenomena should be definitely 
clear in this regime. 

Another interesting issue towards an answer of which, our formalism can 
aim is the study of electrons in a background squeezed state. Squeezed states 
are characterized by negative energy densities, and a connection between 
them and decoherence, dissipation and noise would be very rewarding. 

Finally, the expression (4.10) for the single electron propagator implies a 
spin - momentum self-coupling due to the vacuum fluctuations. The study 
of this kind of couplings is a rather unexplored territory and potentially of 
much interest in the light of quantum optics. 
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A Coefficients in sigle-particle propagator 

We here give expressions for the functions ai, 71 and A in equations (4.17 
- 4.19). 

ai(t,r) = sin^mtsin^mr (A.l) 

47rm 

pJt, r) = /* ds( e-'""^'-'^ smm(s - r) + cosm(2t - s)e^"^(*-^+'-) 

IGnr Jo \ 

-e'""' sin ms - cosm{2t -s + r)e-'"'^^-'^^ (A.2) 
sin TTiT 

7i(t,r) = cos mr + sin m(2f — r) + sin mr) (A. 3) 

IGnmr 

A(t) = At / ds [ ds'(l + cos 2m(s - s') + cos 2m(t - s - s') 
27r"' Jo Jo \ 

+ cos 2m{s — s') cos 2m{s + s') 
+ cos 2m{s — s') cos 2m (t — s — s') 



X7 T7, n (A.4) 

(s - s') + e2 ^ ^ 

Note, that the term A has a cut-off dependence on e, which becomes smaller 
and at late times negligible. It can be absorbed into an early time redefinition 
according to our renormalisation scheme. 

Also, we can see that the term 71 is the one whose amplitude increases 
with time and is to be held responsible for the breakdown of our perturbative 
approximation and possibly as the one giving rise to decoherence. 
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